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ABSTRACT 

Stability  and  convergence  for  a  difference  method  for  quasilinear 
elliptic  boundary  value  problems  are  proved.  Asymptotic  expansions  of 
the  discretization  error,  basic  for  Richardson  extrapolation,  are  estab¬ 
lished.  The  general  theory  of  "discrete  Newton  methods"  and  "iterated 
defect  corrections  via  neighboring  problems"  [6,8]  and  Pereyra’s  deferred 
corrections  [22]  are  used  to  derive  different  high  order  methods.  Some 
special  cases  and  computational  problems  are  pointed  out  and  numerical 
tests  are  included. 
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SIGNIFICANCE  AND  EXPLANATION 

In  this  paper  we  discuss  the  numerical  solution  of  quasilinear 
elliptic  boundary  value  problems  in  a  bounded  region  ft  of  IRn. 

Problems  of  this  type  arise  in  many  physical  and  chemical  applica¬ 
tions,  e.g.  in  chemical  reaction  processes,  in  vortex  motions  of  fluids, 
in  steady  state  heat  conduction,  in  diffusion  processes  and  in  struc¬ 
tural  mechanics. 

By  applying  the  usual  five  point  discretization  in  the  interior  of 
ft  and  high  order  approximation  of  the  boundary  conditions  we  obtain 
stable  discretizations  with  an  error  admitting  an  asymptotic  expansion. 

This  expansion  may  be  employed  with  Richardson  extrapolation  or  defect 
or  deferred  corrections  to  obtain  methods  of  flexible  order  (2,  4  and 
5.5)  . 

Due  to  the  increasing  complexity  of  the  problems,  linear  equations 

* 

often  are  not  adequate,  so  nonlinear  problems  have  to  be  discussed  . 

Error  asymptotic  results  have  been  obtained  before  (Pereyra  [23],  1970) 
for  those  problems  but  under  strong  restrictions  on  the  boundary  of  ft 
which  essentially  are  satisfied  for  rectangular  domains.  For  nonrectangular 
smooth  boundaries  of  ft  only  the  Dirichlet  problem  for  the  Poisson  equation 
was  treated  before  ( Pereyra -Proskurows.ki-Widlund  [24],  1977).  This  paper 
gives  the  first  approach  for  smooth  boundaries  and  nonlinear  problems, 
generalizing  the  results  in  [23,24].  Further,  our  high 

order  methods,  based  on  defect  and  deferred  corrections,  save  a  signifi¬ 
cant  amount  of  computer  time  and  may  be  used  for  adaptive  techniques. 

*Since  only  a  few  special  treatments,  known  in  mathematical 

physics,  apply  to  nonlinear  problems,  one  has  to  treat  most  of  them  by 
numerical  methods. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


HIGH  ORDER  DIFFERENCE  METHODS  FOR  QUASI LI NEAR 
ELLIPTIC  BOUNDARY  VALUE  PROBLEMS  ON  GENERAL  REGIONS 

Klaus  Boomer 

Introduction . 

In  this  paper  we  discuss  the  numerical  solution  of  quasilinear  elliptic  diff»r<i: 
equations  on  a  general  open  bounded  region  ;I  C  jtn  with  vh'.v  .  :  -c:  ; .  ■ 

along  the  boundary  912  (Dirichlet  problem).  Equations  of  this  type  arise  in  several 
practical  problems,  one  of  them  discussed  in  51.  Since  we  are  mainly  interested  in 
high  order  methods  the  situation  has  to  be  smooth  enough. 

We  apply  the  usual  symmetric  five  point  formulas  in  all  recjular  mesh  points  :< 
(all  neighboring  mesh  points  of  x  are  in  Si).  In  irregular  mesh,  print.-  I..  lea-;*, 

one  of  the  mesh  neighbors  lies  outside  of  f!)  we  introduce  the  boundary  con  i  in.-  iv 
using  interpolating  polynomials  of  degree  k  <  •  . 

For  k  =  0  and  1  we  have  the  well  known  first  and  second  order  method  v  tv 
Gerschgorin  (20)  and  Collatz  [14]  for  linear  elliptic  equations.  For  k  y  2  tne 
method  is  due  to  Kreiss  presented  by  Pereyra-Proskurowski-Widlund  (241  for  the  To  is  so 
equation. 

Now,  Wasow  [28]  has  shown  by  some  examples,  that  an  asymptotic  expansion  of  the 
discretization  error  is  not  available  if  the  boundary  values  are  not  reproduced 
exactly  enough.  By  a  modification  of  the  five  point  formula  in  irreqular  mesh  points 
Bramble-Hubbard  [11]  obtain  a  first  term  in  the  asymptotic  expansion. 

Higher  order  asymptotic  expansions  are  only  known  in  two  special  cases:  For 
quasilinear  operators  ft  has  to  be  such  that  there  are  only  regular  mesh  points 
(Pereyra  [23]).  This  regularity  condition  is  satisfied  essentially  only  rcr  recta:.,': 
domains.  For  general  regions  ft  asymptotic  expansions  are  known  only  in  ti.e  iirict.f 
problem  for  the  Poisson  equation  (Pereyra-Proskurowski-Widlund  [24]).  In  this 
paper  we  generalize  the  asymptotic  results  of  Pereyra-Preskurowski-V.'idlur..:  to 
quasilinear  Dirichlet  problems  on  general  regions.  Since  we  are  mainly  interested 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-75-C-0024 ,  the  Deutsche 
Forschungsgemeinschaft  under  Contract  No.  BO  622/1  and  the  University  of  Karlsruhe. 


in  high  order  results,  the  situation  has  to  be  smooth  enough.  Our  results,  a;  .  r 

Wasow’s  (28),  are  obtained  here  by  combining  stability  and  local  discretization  errors. 

Ir.  [10 J  we  use  special  techniques  to  rive  far  linear  elliptic  :  ~i  ■■  ■ 

considerably  better  results.  Especially  the  conditions  which  we  have  to  impose 
the  coefficients  are  much  less  stringent  than  here  and  the  error  asymptotic  for  ■ 
is  nearly  as  good  as  here  for  k-\ . 

In  4;  we  formulate  the  iiscreti.tatio::  schem,.-  an:  describe  -no  equation,; 
obtained,  51  is  devoted  to  a  stability  proof  fir  our  method  for  k<6  using  some  o: 
lemmas  or  treiss,  given  in  [24].  Since  we  have  to  impose  in  ;  •  relatively  strong 
conditions  on  the  differential  operator  obtaining  only  stability  with  respect  to 
the  >.,,-r:orm,  we  give  another  stability  proof  in  '4  ,;u-r  mum  weaker  --n  ft! 
for  the  Za-n orm,  but  only  for  k<4.  Tor  k>0  we  can’t  prove  •  ,f  i.ity. 

The  convergence  of  the  basic  second  order  method  and  the  existence  of  an  asymptotic 
expansion  of  the  (global)  discretization  error  are  proved  in  15  for  k  <  6.  Two  special 
cases  which  allow  an  extended  ••asymptotic  expansion  are  discussed  in  §6,  one  of  them 

f. 

treated  already  in  connection  with  iterated  deferred  corrections  by  Pereyra  [23) . 

§7  uses  the  asymptotic  results  in  §  § 5  and  6  to  obtain  high  order  methods  via 
Richardson  extrapolation.  Further  it  combines  the  asymptotic  expansion  with  the  general 
theory  of  '(discrete  Newton  methods"  [8]  again  providing  high  order  methods.  In  addition 

we  indicate  high  order  methods  based  on  iterated  defect  corrections  via  neighbouring 
problems  [29,  30,  27,  15,  16,  7,  8]  and  on  iterated  deferred  corrections  formalized  by 
Pereyra  1221,  The  mutual  advantages  and  disadvantages  of  defect  corrections  and 
Richardson  extrapolation  are  discussed.  Especially,  when  ail  mesh  points  are  regular 
the  asymptotic  expansion  nay  be  extended,  depending  essentially  on  the  smoothness 
of  the  solution.  Ln  that  case  Richardson  extrapolation  is  not  able  to  compete 
with  the  other  methods,  especially  with  the  discrete  Newton  methods. 

The  discussion  of  several  problems  arising  in  the  computation  of  the  different 
methods  are  the  topic  of  ‘,8.  The  applicability  of  fast  Laplace  solvers  (see 
Pereyra -Proskurowski-Widlund  [24]  and  of  Newton  methods  for  the  original  nonlinear 
discrete  problem  are  discussed. 
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1 .  A  quasilinear  elliptic  boundary  value  problem  in  chemical  physics 

In  a  chemical  reaction  let  T  be  the  temperature,  k  the  thermal  conduct iv: •  , 

Q  the  monomolecular  heat  of  reaction  and  V  the  velocity  of  the  reaction.  Then  tc 
equilibrium  of  the  heat  transfer  between  the  heat  produced  by  the  chemical  reaction, 
and  the  heat  conducted  away  is  characterise  :  in  ..  e  by  the  :  > .  •  [  ; 

p.  Is 01: 


(1.1) 


kAT^-.V  it  n  (A  :=  I  ~  )• 

i=l  ?x 


Now,  V  and  T  are  related  by  the  Arrhenius  relation,  which  reads,  in  the  simpl-. 
case  we  have  assumed  here,  as 


(1.2) 


V  =  cv  exp(-  —  )  , 


where  c  is  the  concentration,  v  a  scaling  factor  which  may  be  positive  or  negative 
E  the  energy  of  activation,  and  R  the  universal  gas  constant.  So  we  obtain  from 
(1.1)  and  (1.2)  the  equation 


(1.3) 


kAT  =  -  cvQ  exp(-  —  )  in  Cl 
RT 


with  suitable  boundary  conditions  on  3C1 


For  constant  Q  and  k,  and  T  close  to  TQ,  (1.3)  may  be  transformed  into 


(1.4) 


A\Ji  +  u  exp  tp  =  0 


.  C-V.Q.E  ,  E  E 

with  u  ==  - 5 —  •  exp(  )  ,  $=  - s 

kRT  0  RT  ‘ 

0  o 


(T-T  ). 


It  is  possible  to  give  the  general  solution  for  (1.4)  in  the  two  dimensional  case 
(n  =  2  in  (1.1)),  see  Ames  [2],  p.  182.  The  approach  given  there  may  be  generalized 
to  the  equation  (1. 3) .  To  obtain  a  solution  satisfying  the  boundary  conditions  in 

(1.3)  one  has  still  to  solve  the  difficult  problem  to  choose  three  suitable  fur  :tio'.: 
F,f,g  such  that  i|>,  defined  by  exp  ip=F(f,g), satisfies  the  boundary  conditions!  tee  [.' 

For  n  >_  3  or  for  anisotropic  heat  conduction  this  approach  is  no  longer  possibl 
and  we  have  to  find  other  methods.  Here  we  are  concerned  with  the  corresrcr. : ' 

difference  analoga  to  (1.3)  and  their  generalizations  and  with  improvements  of  the 


approximate  solutions. 
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2.  The  basic  difference  scheme  for  elliotic  problems. 


Let  a^,  i  =  l,...,n,  f  and  g  be  continuous  real  valued  functions,  do f i r.c  1 


a:,  open,  bounced  and  connected  f.  cF,  . 


resoective-v 


'  ..hen  we  cefme  an  elliptic  differential  operator 


(  a  .  i  •  i  ■  ■■  >  7y 


[d  :=  {ye  C IP.)  <"i  C(.Q)|  <y(x),7y(x))  e  B-  -*  CC)  >  C  O  ) 


(2.1)  F 


y  Fy 


I  a.  <•) yx  x  (-)  +  f ( • , y ( • ) , Vy ( - ) 


I  i=l  i 

[>■(•)-  g(-)  in 


{a.  (•)  ■»  a  >  0  . 
1 - 


Now  let  Fy  =  0  have  a  unique  solution  z  in  D  with 

K  ••=  (y  e  d|  lly  -  z ||  L  (S)  i  o)  , 


2  is  unique  solution  in  D  of 

z 


-  £  a  (-)z  (•)  +  f  (•  ,z  ( • )  ,  7z  ( • )  )  =  0  in 


0  =  F2  =  i=l  1  Vi 


l  |z(d  -  g(d  =0  in  3£  J 

Conditions  for  unique  solutions  are  given,  e.g.,  in  Bers  [3], 

To  define  a  discretization  we  introduce  a  grid  !\  and  the  grid  liner  Jv  _  a: 


:  =  {x  e  ®n|x  = 

(X. 

h,  n 

*h,n 

=  {x  e  JRn  |x  = 

(x. 

I  n  T 

IG,  :  =  {xe]R  jx  =  (x, , ,x  )  ,x  e  E,x.  =n.h,n.  f  IZ  for  i  ?  v  :  ' 

v  h,n  ‘  1  n  v  l  li 

where  we  have, for  simplicity,  chosen  an  equal  stepsize  h  in  all  directions.  The 

following  discretization  is  a  generalization  of  methods  given  by  Gerschgorin  [2  1  and 

Collatz  [14]  for  linear  elliptic  second  order  equations  (k  =  0  and  k  =  1)  and  by 

Kreiss  and  Pereyra-Proskurowski-Widlund  [ful  (k.  <  6)  for  Poisson's  equation. 

With  the  open  set  (2  we  introduce  the  mesh  points  in  .1  and  distinguish  between 

the  set  of  regular  mesh  points  (slightly  different  from  the  Introduction) 

(1  :=  fx  e  Sllx  +  he,  e  n,  i  =  l,...,n)  n  r 

h  —  l  h,n 
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(here  e  are  the  unit-vectors  in  the  direction  of  the  positive  i-th  coordinate  axis) 
i 

and  the  set  of  irregular  (mesh)  points 

:  -  i  x  f  -  i  3  i  :  x  +  he  /  or  x  -  he.  /  ^  r 

h,  1  i  l  n ,  n 

In  regular  points  we  use  the  standard  centered  difference  approximations  for  the 

derivatives  to  obtain  ...  .  *.  r;vt  ;  .  1:.*  :I  •  i  i  *?  !'  the  fort 


2  r  . 

h  (vF)n(x):=  -  )  a.  (x)  \n,  (x  +  he.)  -  2n  (x)  +  n  (x  -  he . )  r 
hh  -  h  l  h  h  l 

i=l 


(2.4)  \ 


_  n.  (x  +  he  )  -  n.  (x  -  he  ) 

.2  ,  ,  h  .  h 

h  f  ( x ,  "i  ( x  > ,  - — - 

n  2n 


1,  (x  +  he  )  -  i,  (x  -  ho  )  ] 
h  n  h  n 


2h 


In  this  case  we 


for  X  t  . 

h 

If  x  is  an  irregular  mesh  point,  at  least  one  of  the  x  +  he^  1 

have  to  replace  in  (2.4)  every  n.(x  +  he.)  with  x  +  he.  f  .;  by  a  provisional  value 

h  —  1  —  l 

obtained  by  polynomial  extrapolation.  Let,  for  that  purpose,  x  e  D.  t  ... 

h,  i 

but  x  -  he.,..., x  -  (k  -  l)he.  e  il,  .  Further  let  xf,  x  <  xf  <  x  +  he.,  be  the 
i  i  h  i  l  -  l 

unique  intersection  of  and  the  line  segment  [x,  x  +  hej.  For  h  small  and 


smooth  enough,  these  conditions  will  be  satisfied. 


s .  h 

i 


Now  we  define  the  following  approximations  for  y  (resp.  y  ):  Compute  an  inter- 

i  i  xi 

polatinq  polynomial  P  of  degree  k,  defined  by  P,  (x  -  (v  -  l)he.)  =  y  , 
k  k  l  l~v 

2 

v  =  l,...,k,  P,  (x?)  =  g(xf),  and  replace  inhy  ssy-2y+v  resp 
k  l  l  x.x.  'i  yo  y-l 

li 

2hyx  »  yx  -  y_1  the  value  yx  by  :  =  Pk(x  +  he.).  If  the  distance  from  x*  to 
x  +  he^  is  s^h,  with  0  <_  <  1,  we  find  (see  Pereyra-Proskurowski-Widlund  124]) 
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Again  by,  e.g.  Stetter  [26],  we  need  for  the  proof  cf  the  stahilitv  of  r f. •: 
(<P^F)  only  to  discuss  the  stability  of  ((p.  )  ft;-  '  t.  -A.  c  sr.aii  er.cn,-:. 

(2.h)  and  (2.6)  imply  that  ,  for  f  6  CJ'+'1(S)  arc  v  €  with  esset.t i.-iliv  : 

derivatives,  «^F’(y)  =  (oyV.iy)  +  0(h‘).  That  means,  we  may  restrict  cur 
proof  to  the  discussions  of  ip^FMc^):  So  we  formulate  the  discretisation  ip 
affine  operator  of  the  form  (2.1).  Kith  continuous  real  valued  functions 
ai^  ’  ^i^  ’  i  =  1»  •  •  •  ,n,  c(  •  )  and  d(  •  ),  defined  on  C.,  a.  ( •  l  ■>  a  >  :  a:  ; 

g( •  )  defined  on  30  we  have 


In  irregular  points  those  nh<x  +  he.)  in  (2. 10)  have  to  be  replaced  by  expressions 
corresponding  to  (2.5)  for  which  x  +  he^  ^  For  the  stability  proof  in  ?3  it  is 
only  important  that  <0^,  FQ  :=  Fd|d=0>  may  be  separated  into  the  sum  of  dis¬ 
cretizations  for  ordinary  second  order  equations.  So  we  have  for  regular  points: 
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(2.11) 


h  (V,  F  )n  (x) 
nu  n 


V  i r  2 

l  a.  (x)  2  +  h 


c  (x) 


1=1 


-rh(x) 


L  a, 

v=l  ' 


,  b . (x) 
-1  ♦  *  1  - 


2  a .  (x) 


■  n.  (x  +  he.  )  +  i  -  1  - 
h  l 


,  b  Ax)- 
h  i 

2  a  Ax)'  "h 


(x  -  he  t 


Again,  for  irregular  points  the  r»  (x  +  he.)  in  (2.11)  with  x  -*■ 

n  ~  i 

replaced  by  expressions  corresponding  to  (2.5).  Further,  we  add 


points  tc  obtain  a  oyster 


Br.uugh  to  study  the  properties  of  the  :r.atr:x  A  ;e:: 


(2.12) 


h  (v'hF0)nh(x) 


for  x  e  '  and  for  x  c 
h 


This  matrix  A  may  be  written  in  the  form  (see  Pereyra-Proskurowski-KidU 


(2.13) 


>  P.A.P. 

.  L,  111 
1  =  1 


with  suitable  permutation  matrices  P  .  The  matrices  A  are  obtained  by  col  lectin- 


only  those  contributions  in  (^,F„)n,  (x),  which  are  multiplied  by  a  fixed 

h  u  h 


a  .  ( x )  ( :■ 

l 


(2.11)).  Therefore  these  A^  are  direct  sums  of  matrices  of  the  form  (  . - 
we  use  the  abbreviations  to  be  introduced  in  C.-i) 

(2.14)  B:= 


ai  +  -  a2  +  a3  + 


/  2  a.  ,  ^  . 

2+h  q  +— rP  »“P  +  ~ rP  »  — TV 
Ho  <x'*o  *o  a'  o  a*  o 
o  o  o 


1  ! 


"Pi 


,  2+h  q1  ,  -  ,0  ,  0  , . 


+  ,  ,2  -  _ 

-p  ,  2+h  q9,-p9  ,  0 


“k  +  n 

’a^o’  0  • 

o 


0  ,  C  ,  . 

.,0,0 . 


4-1 


,  0  ,. 

o 


-pc. 


?c+r^t’ 
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[x  ,x^]  C  J  "  g,  in  the  direction  of  the  coordinate  vector  ■ 
i  1  h,n 

see  Figure  2.  For  simplicity  we  have  used  the  following  abbreviations 


b.  (•) 

a(  •  )  :=  a .  ( ■  )  |  ,P<->  :=  ~~T  i  . 

1  1 |x . , x*  I  ai  !  lx , , xf ) 

li  ii 


.  q<->:  = 


i  '  •  1  . 

v=l  ’  L  .  ,  •  *  ' 


where  av,i  corre-?c:i  t0  ’ne  -rs: 

i *.tersect  L  :-:t  or  c>.~.  a:.  :  {x-te .  '  t >0>f  lef  t  o:  x*.  further 

let  a  :=  a<>: :  ),  p,  :  =q<  1 ,  where  are  :o.i:.e:  ry  ."iyur-  . 


+  h 

V  =  *  +  7  'v  ’  V  =  1 


Figure  2 

In  our  discussion  we  have  confined  the  coefficients  a. (•)  to  be  functions  of  the 

1 

independent  variable  x  €  ft  alone.  Yet,  we  could  generalize  the  whole  discussion  in 
§§2  ff.  to  the  case 

n 

-  £  a.(*,z(*),z  (*))z  (•)  +  f  (•  ,z(* )  ,Vz(* ) )  in 

(2.16)  Fz  :=  i=l  1  Xi  Xi  i 

z(  )  -  £(•)  on  3.Q 

with  a^f-  ,z(*)  ,z  (.))  >  a  >  0. 
i 
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Then  we  would  have  the  derivative 


I  '  a  .  ( - ,  y  ( • ) ,  v  ( -  >  >  u  <•)  <■  a10’1,0’  (-,v(.),v  ( -  >  >  u  ( - 1 

1  X  XX  1  X 

1=1  111  1 


(2.17)  F* 


*  (y )  u  =  ■!  +  a .  '  )f  ( • ,  v  ( • ) ,  v  ( •  )  u  •  ♦  f 
l  *  x  x 


1  5)  l-,y(-l,vv(.))u(-l  \ 


I*  i  f10’---’1 . ■' 


(*,v(»),  v(.))u  (•)  in  u(-)  on 

x 

i 


which  is  of  the  same  structure,  but  more  cor;  licated  than  ' 


In  53  we  qivo  a  stability  proof  for  (2.12)  and  k  (  by  showing  that 
!! A II 2  _  Ch  ^  for  the  matrix  A  defined  in  (2.12).  So  all  results  based  on  5  3  are 
results  for  the  Euclidean  norm.  The  conditions  which  we  have  to  impose  are  rather 
stringent.  Therefore  we  qivo  another,  much  easier,  stability  proof,  valid  only  for 

k  <_  4  showir.q  that  || A 1 1 ^ Ch  So,  using  this  ?4  with  its  less  stringent  conditions 

•  •  •  •  «••••• 

sup-norm  results  are  obtained.  Since  k  essentially  -•  ie:’ 

in  the  asymptotic  expans  ior:  u:.  i  ter.r  •  ire  w  u  ti.w:. :  . 
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3.  Stability 

By  stability  we  mean  the  uniform  boundedness  of  the  inverst  o:  •  r.i*  _ :  * 

the  difference  operator  in  a  r.eirhh turr.c  .  :  :  t:.*-  a  *  *:.■ 

To  prove  this  property  for  the  method,  defined  in  (2.4>-(2.7),  it  is  <  nou  t 
stability  of  the  discretization  for  the  linear  problem  F_  in  r 

F  y  :=  F 1  (z) y  a:.:  f  from  tee  tetter  [.  ]  ,r .  *:  i.  ,  .  :v  ■ 

prove  that  the  matrix  A,  defined  by  (2.12),  has  the  property  A  ^  r. 

For  that  purpose  we  use  and  generalize  some  of  the  Lemmas,  due  to  Kreiss,  :iv» 
Pereyra-Proskurowski*-Widlund  [24].  The  norms  in  this  paragraph  are  <  x  ;  :-:v*  . 
Euclidean  vector  norm  and  the  spectral  matrix  norm 
Lemma  3.1  ( 24 )  :  Let  the  symmetric  part  of  a  matrix  A  sat i sf y 

(A  +  AT)/2  >  51,  5  >  0  . 

Then  A  is  regular  and  ||a  ^  [I  <  6  L 

5  T 

Lemma  3.2  [241;  Let  A  =  £  P^A^P^  with  permutation  matrices  F^.  lf_ 

•  •  ••  ••  .  •  •  •• 

(A.  +  AT)/2  >  «I,  f  >  0  ,  for  ail  i. 

i  x  - 


Then 


Lemma  3.3  [ 24 ] : 

(3.1) 

Then 


(A  +  A  )/2  >  nS I  . 

Let  the  matrix  A.  be  the  direct  sum  of  certain  matrices  B.  .  . 

- - - — — — -  j_ - - -  n 


(B. .  +  B. ,)/2  >61,  6  >  0,  for  al]  j  , 

13  13  -  - - 


If 


(A.  +  AT)/2  >  6 1  . 
11  — 


We  have  seen  in  §2  that  the  matrix  A  i:.  (2.1!)  i?  the  sum  c:  matrices  .. 

the  A.  are  direct  sums  of  matrices  of  the  form  (2.1^)  ar.d  r.  are  oermuta t 1 matr:  ■ 
i  i  c 

we  have  assumed  h  so  small  that  x  e  ^  p  implies,  e.g.,  x  +  hei €  C  and 

x  -  (v  -  l)he^  €  for  v  =  l,...,2k  -  1.  Especially  we  have  therefore  in  (2.14' 

that  l  >  2k  -  1 .  With  the  functions  a,  p,  q,  3rd  a  ir.tr  ::  \  .  . 

v  v  v 

with  indices  changed  for  convenience,  vie  may  prove  as  ir.  (.  ■»]  the  foilovsir.r 
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I 


Lemma  3.4:  Let  the  matrix  B  in^  (2.14)  be  split  into  two  matrices  B.  ,  P.  ,  whvi>.  . 
with  p*  and  p'  in  (2.15), 

(3.2)  B  := 


Let  be  obtained  from  by  inverting  the  order  of  rows  and  columns  and  changing 


m  to  m1  and  a ^  to  ah  Further  let 


If 


m  +  m’  =  l  '+  l,*  m  ■>  k',  *m,*>  k  .  *  • 


(B1  +  b’J’j/2  >  «I  and  (B2  +  B^/2  >  <51  , 


then 


(B  +  BT)/2  >  61  . 


Before  we  are  able  to  prove  the  next  Lemma  we  have  to  do  some  preparations.  We  obtain 


the  symmetric  part  S  of 


(3.3)  S  :=  ( B1  +  B^)/2  = 


0 

0 

Smm-k+l 

s  , 
mm-i 

s 

mm 
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where 


(3.4) 


s  =  a  (1  +  —  q  ) ,  s  ~  a  (2  +  h  q  )  ,  v=2,...,m-l, 

11  1  2  1  vv  v  v 


mm  m 


v+lv 


1 2  +  h  a  +  (1  *  -  p  )  —  » 

V  m  2  m  a0J 

a  ,  +  a  , 
v+1  v  h 

-  +  —  { a  p  -  a  p  ) 

2  4  v+lpv+l  v  v 


s  =  +a  (1 


mv  m  2  m  2a 


P  > 


m-y+1 


v  =  2 ,  .  .  .  ,m  -  2  y 


v  -  m  -  k  +  1  , . . . ,m  -  1  , 


mm-1 


a  +  a  _  .  .  a_ 

m  m-1  n  ,  .  , ,  h  ,  2 

-  +  —  (a  p  -  a  p  )  -  a  (1  +  —  p  )  - 

2  4  nr  m  m-rm-1  m  2  Fm  2a  J 


s  =0  elsewhere 
vp 


Let  us  now  assume  that 


(3.5) 


r  3+a,  +  *  2+a  +  *  1+a  +  * 

1  a*C  (x . / x . ] ,  P€C  [x. ,x. ] ,  q  e  C  fx. ,x. ] ,  0  <  a  <  1  , 

li  ii  ii 

+  * 

a(-)  >_  a  >  0  in  ,  x  ^  ]  .  There  exists 
-  +  +  »  * 

x  *  [x^  +  (k  +  1  *■  s^h,  x^  -  (k*  +  1*  -  s^l+i]  with  •  • 

a'  (x)  =  0,  a'  (•)  >_  0  in  (x*\x)  and  a'  (■)  £  0  in  [x.x^J  . 

a'2]  +  * 

2aq  -  a"  t>. sgn  a')(ap>  +  —  <•)  >  c+  >  0  in  [x^x..]  , 


'  21  * 

2aq+(sgn  a')  (ap)'  -  -r—  ( • )  >0  in  [x.,x.]  .where  we  define 

Za  )  11 

{2aq  -|(ap)'|}(*)  >c+>0  in  an  interval  of  length  >  h  with  the  midpoint  x. 


sgn  a' 


1  in  (x^.x)  for  a  ^  0 

/V  * 

-1  in  (x.x^]  for  a'  <  0 


Here  Cm+ot  Ixt,xf]  means  that  all  derivatives  up  to  the  order  m  are  (uniformly  )  Holder- 


continuous  in  tx^x^]  with  exponent  a,  0  <  a  <  1.  For  |y  •  - a"  j  ( ■ )  <  c+ 

the  forelast  inequality  ( >c  + )  in  (3.5)  includes  the  last  (>0).Ke  will  come  back  to  the 
a  '<•)  <  0  in  [x*,xj,  a'(-)  ^0  in  [x.x*]  in  Remark  3.7. 

Now  we  choose  the  numbering  in  (3.2)  such  that  ]x-x  !<  - 

1 1  _  2  • 
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(3.5) 


where  the  Z  and  the  r  are  defined  in  the  usual  way,  so,  e . q . , 
uv  v 


ri  :=  sn' 


(  3.7) 


v+1 


V+1  V 


V+1  V 

*  s 


V+1  V 


?  ,  /  r  , 

V  +  1  V  V 


v+1 v+1 


v  =  1  ,  .  . . ,m  -  2 
V  =  1 , . . . ,m  -  2 


To  give  the  Z  and  r  we  would  have  to  introduce  the  whole  algorithm.  Since  we 

•  •  •  •  •  MV  •  •  HI  .  «  •  •  •  •  •  •  •  *  • 


need  only  the  properties,  given  in  (3.13)  for  the 


!  and  in  (3.34)  for  r 
mv  m 


we  do  not  give  the  explicit  results  corresponding  to  (  .). 


Using  the  symmetry  of  S  and  the  properties  of  the  l  one  stral.-ht  torew . 

verifies  v+lv 


(3.8) 


S  =  L 


So  we  will  study  the  matrices 


(3.9) 


and  D 


r 


1 


o 


r 

m 


we  again  have 


(3.16) 


0  for  k  =  h  small  . 


S  >_  61  with 

<S  =  Cd  .  h2/(di am(B) )  2 
nun 


Proof:  Since  we  need  estimates  of  the  form  (3.17)  we  give  lower  estimates  for 

T  T  T  ^  1 1  T  M  2 

x  Sx  =  x  LDL  x  >  min{r  }  •  L  x  ,  . 

—  ,  v  "2 

v=l 

1 

Now  by  (3.5)  and  with  v  <_  m  <_  0  (^ )  we  find  by  (3.13)  +  0(h),  so 

>  a ^/2  >  a/2 ,  v  =  1,2, ... ,m  -  1,  and  with  (3.16) 

(3.18)  xTSx  >_  min{a,c+  (k)  }/2  •  ||lTx||2  . 

— T  _rr  (  f 

We  try  to  find  an  upper  bound  for  j|  L  y  ||  0  since  1  y  f  j0  <  cl!  y'|  ,,c  €  f.  ,  i-;pl:«- 
II  LTx ||  2  >  c-1 1|  x  ||  2  .  With  ym  €  F,  yT  =  (yT,ym>  and  £‘,  Ln  in  (  3 . Q ) ,  one 


verifies,  since  L  is  regular, 


yT  L'X  =  (<yT  -  VKI'V  - 


simply  by  multiplying  the  equation  from  right  with  L.  Since  the  “v/a0  are  of  the 

s 

6  _(k)  •  -  and  since  the  vector  £  in  (3.13)  includes  only  a  /a  ,  v  =  2,...,k, 

V  /  U  V  —  S  VO 

has  a  uniformly  bounded  norm  in  0  <  s  £  1.  So 


lL  Ty II 2  —  II Lii tl 2 ( ^ yT 


<lliTll,  ♦  lyJIWI,)2  a  I y  | 2  i  c(||l"| ||2  ♦  i)||y| 


Now,  since 


I  Lix  II 2  =  ^smallest  eigenvalue  of  L^L^)  1  , 


we  are  going  to  estimate  this  eigenvalue.  The  definition  of  L  in  (3.6)  and  (3.9) 


"2]  1  *  '  21  *  32 


T 

L  L 
11  11 


?  ]  -f  f 

v .  - 1  vv-1  V+  I  V 


\  \ 
\ 


n-lm-2  ''m-  lm- 2  j 


We  use  the  well-known  theorem  of  Gerschqorin  [18]  to  : : 


.  -1  ,,  ;  j'  ’  '‘'21!'  1  *  'm-lm-2  "  '  'm-lm-2  )  ,2 

|L,  ,|l-  1  mm2  f  ;•  h 

U  2  U  +  t2  ,  -  i  v  =  2, - m  -  2 

v  vv-1  vv-1  v+lv 


with  c ^  €  3R+.  We  find,  with  'x.-x'  <  — ,  a"(x)  <  C  ar.i  *-v+<v=  -l+2(h) 

/l-j  iL,  !  >  —  {2aq-(ao)’ }.  for  x.  <  &, 

1  d1  -  *  i  -  - 

.  2 

(3.20)*  1- !  ^2i !  -  {2aq-(ap),-aM  }.  for  x;  >  ft, 

•  1  •  .  •  .  .  • 

i +l2  ,  -  -|H  .  =  i+Cfh ) . 

\  m-lm-2  1  m-im-2 

Ignoring,  for  a  moment,  the  vC( h" )-terns  in  (3.12)  we  find 


1  +  1  ,  -  K  , 

vv-1  vv-1 


2  ,  -  (1  +  f  ,)•((>  ,  -  f  ,  ) 

v  +  lv1  vv-1  v+lv  vv-1 


-  T  {Jq  '  >f>  - 


,3  a  v-2  .2 

+  —  l  (2  (ap) 1  +2 - 4aq)  *  o(h3)  . 

2  a2  ,  u-1  3  " 

v-1 

Since  V  might  be  I  we  cannot  neqlect  the  h3  •  T.  -  term.  Combining  (3.20)  and 

(3.5)  we  find  (3.19)  satisfied.  We  have  ignored  the  vO(h3)  -  terms  in  f  resp. 

v+lv 

(v  -  l)0(h3)  in  That  we  may  do  so  is  shown  by  an  elementary  straightforward, 

but  very  lengthy  computation,  essentially  by  showing  that  for  these  vO(h3)-terms 
2  (v  -  l)0(h3)  -  vO  (h 3 )  -  (V  -  l)0(h3)  =  0  (h3)  . 
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As  a  consequence  of  (  3.10  ),<  :  .  .i  1  )  we  nave  prove:  (  ■ . '  *) 
(  3 . 1 S )  Lemma  .i' .  3  .  □ 


in.:  mere r ore  w : t r.  t. 


Remark  t.r;  The  conditions  in  (  are  only  cr.e  possibility  t;  enforce  (  j  .  i 
only  purpose  for  a'  <  0  and  iaq-(a*  ^  '  :  j  to  ensure ,  that  the  h 

term  in  (3.2i)mmt  not  dominate  the  im-terr  in  (  .'r.e  Ti;ht  choose  oth 

limiting,  properties.  We  get  then,  with 


a  :=  (a]))  -  (a{  )  .  ♦  j  ( —  -  4aq)dx 

v~  1  v- 1  0  a 


condition.",  o t  tne  lorm 


f  a*2)  v“2  ..2 

a  *2aq  -  a"  -  (at.)’  ♦  — — j  +  a’  Y  (2  (a t>)  ’  +  J -  -  4ao) 

^  ‘  a  •  v  £*i  a  i 

v  u  =  1 


wher*  i-'v-  h  or 


i  U  -  y  *  -ir1)  *  m  Vi  i  cf 


U  -  fe-)’ 


isEii)  4  _i_  q  (a.  .  , ,  c 

a  *  2  9v-J  v-1  8  -  «■ 

v  a  , 
v-l 


or  MTiatioi:  c:  *thc*:.e  ;m:j  r  ions  w?  *.?.  1 1  :.-ur :  ce  .  part:  li  iiVcv 

oqmti  1  r*ven  more  cor.f  1 1 :  tied  to  check  :r.  con  problems  than  (  -  ,r  ) 

w*  coniine  our  further  ii  ;  :!-r.  to  i  ;  f  in  a  cpcciil  ..ice  (  • . 1  )  rd.oul 

satisfied  one  still  can  try  to  verify  the  -orrespond i ng  conditions  based  on  * 
inequalities  given  in  (.<..*).  □ 


Remark  3.7:  We  have  confined  our  discussion  to  a*  (•)  >.0  in  [x.,x],  a’  (*)  <  0  in 

„  * 

[x,x^] .  If  the  distribution  of  signs  changes,  we  have  to  impose 


_  -3+Cl  -  +  *«•  c  r  +  *  ■>  -  ,1+Ci  r  +  *  1  -  ,  _  ,  . 

a  €  C  ixi,xiJ,p  €  C  bti,xiJ,q  €  ^  [x^x^J,  ^  <  .x  <  . , 

+  * 

a(  -  )  >  a  >  0  ir.  bt.,x.  ]. 

^  +  +  *  * 

There  exists  x  €  [x.  +(k+l-s .  )h,x. -(k+i-s .  )h]  with 
1  i  i  i 

a'(x)  =  0,  a'(a)  <  0  in  [xt,x],a'(’)  >  0  in  [x,x^]  , 

( 3 . 23  )j  {2aq- 1  (ap)'  |  }( •  )<C_<  0  in  a  neighbourho:  d  of  x, 

a  t  ^  +  * 

{2aq-a"-(sgna,)(ap) '  +  -  }(•)  <  C  <  G  in  fet.,x.], 

a  -  -  ii 

afZ  +  *■» 

teaq-Csgna'Xap)* - —}<  • )  <  0  in  tx.,x,J  , 

a  -  ii 

{/».  * 

1  in  (x,x.]  , 

+  1 

-1  in  k^,x). 

Then  we  find  S  <  -  61,  6  >  0  and  Lemmas  3. 1-3. 4  remain  valid,  if  5,  >  is  replace 

by  -6,  <  and  the  main  result  |j  A  1  j|  <6  1  stays  unchanged.  So  we  find  again  the 

—  1  —2 

stability  result  ||  A  ||  <  Ch  ,  C  >  0.  □ 

Remark  3.8:  If  we  want  to  prove  the  result  corresponding  to  Lemma  2.5  for  1 ,  the 

order  of  the  indices  in  (3.2)ff  has  to  be  inverted.  That  may  be  simulated  by  the 

transformation  of  x  into  -x,  and  ir.  (3.20),  (3.21)  a'  and  (ap)'  change  into 

. .  •  •*  - 

-a'  and  -(ap)',  the  rest  stays  unchanged.  So  conditions  (3.5)  imply  Lemma  3,5 

for  S  defined  by  and  by  .  □ 

Finally  we  want  to  point  out  a  very  important  special  case: 

Remark  3.9:  If  a  =  const.  Lemma  3.5  remains  valid,  if  (3.5)  is  replaced  by 


(3.24) 


p  £  C2+a  kt.xj  ,  q  €  C1+alxt,!<i],  a  >  0,0  <  a  <  1 
(2q-|p'|)(>)  >  c+>  0  in  kt.x^l  . 


In  (3.24)  one  may  replace  2q-|p'|  by  one  of  the  conditions 
2q  +  p'  for  x  >  x  and  2q-p'  for  x  <  ft 

respectively 

2q  -  p'  for  x  >  ft  and  2q+p'  for  x  <  ft 

4  +  *  * 

where  8  is  an  arbitrary  point  in  t(k+l-s^)h,x^  -*(k+l-s^)hj.  □ 
Combining  Lemmas  3. 1-3. 5  we  obtain 
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Theorem  3.10:  Let  a ,  p ,  q  introduced  in  (2.15)  satisfy  one  of  the  conditions 
(3.5),  (3.24),  (3.23)  or  (3.22)  ,let  k  <  6  and  h  be  small  enough,  and  let  A  be  the 


matrix  defined  in  (2.12)  -  (2.15).  Then  there  exist  constants  C.  >  0,  independent 

-  -  k  - £ - - 


A-1 1|  2  <  Ck(  diamft)2  •  h2. 


h,  such  that 
(3.25) 

If  we  apply  Theorem  3.10  to  a  linear  elliptic  equation  of  the  form  (2.9!  we  obtain  < 
stability  result  in  the  sense  that  the  inverse  operator  to  *  defined  in  (2.:r 

is  uniformly  bounded  for  h  •*  0.  For  the  nonlinear  discretization  (2-4)- (2.7)  the 
stability  is  understood  in  the  above  sense  for  the  operator  F'(y),  where  y  is 
close  to  the  exact  solution  z  of  (2.2)  (see  Stetter  [26]).  So  this  "nonlinear" 
stability  property  strongly  depends  on  small  enough  ||y  -  z||.  We  only  formulate  a 
stability  result  in  this  nonlinear  setting  for  problem  (2.2).  We  will  treat  the  more 
complicated  problem  (2.16)  with  the  much  less  restrictive  conditions  in  §4. 

Theorem  3.11:  Let  k  <_  6,  h  be  small  enough,  and  for  small  ||y  -  z]|  and  i:.  ( .. . ; 


(3.26) 


q(-)  s-  f 


(0,1,0 . 0) 


1  =  i,2,.. 

n 

( • ,  y  ( * )  ,  vy  <  • )  )  /  l  a  (•)  e  c:  <P.),  0  <  -< 

v=l 


4  *  — 

Further,  let  for  each  maximal  intersecting  interval  [x^,x^]  C  n  g  J 

"t*  ^  ^ 

x.,x.  €  3ft)  in  the  direction  of  the  co-ordinate  vector  e.  exist  >:,.c 

i  *  i  - - - —  i - 

+  +  *  * 

[x.  +(k+l-s.)h,  x.-(k+l-s. )h]  such  that 


a.(x.)  =  0,  a.(m)  >0  in  fx.,x.},  a.(-)  <  0  in  fx.,x.J 
11  i  —  —  ii  i  —  ii* 

{2a.q-|(aip.)' |}(. )  >  c+  >  0  in  an  interval  of  length  >  h  with 

^  middlepoint  x, 

a .' 2  t 

{2a. q  -  aV  +  (sgna.')(a.p.  )'  +  ——}(•)  >  c  >  0  in  [xT,x.] 

i  i  ill  a .  _  T  l  i 

a!2 

{2a. q  +(sgna|)(a.p.  )'-  -r^-H  •  )  >  0  in  [xT,x.]. 

i  i  iri  2a^  -  ii 

Then  the  discretization,  defined  in  (2.4)-(2.7),  is  stable  in  the  sense  discussec 
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above.  Instead  cf  (  3 . 1  ^  )  me  coni.:  ief  ir.e  the  rc-rreepcr : ir.g  c:. c i t i crs  ba ? v  : 
or  more  complicated  conditions  base!  or.  (  ?  ,  I 1  ) . 

Vfe  want  to  formulate  t'r.e  corresponding  recnl .  imr  constant  a.. 

Theorer  !•.!,;  bet  V  <  0 ,  e  sr.a  11  ere  ugh ,  a.  -  c-T.st  an.:  rcr  or, all  er-zecr. 
v  -  z”  and  sr.ee th  ereezr  v 


P .  <  - )  :  =  f  ‘  ‘ .  ,y  ( • )  ,  Vy  ( • )  )  f  C2  +  1(  ) 

1  i  -*-2 


(3.28) 


q(-)  :=  f<0,1,0'  01  (-.y(-),Ty(-n  f  c1+',(.:i. 


n  •  ic] 


2q  ( * )  -  7  p >_  c #  •*  0  in 

i  =  l 


Then  the  discretization,  defined  in  (2.4) -  (2.7)  is  stabl e . 

Proof ;  Theorem  3.11  is  an  immediate  consequence  of  Theorem  3.10.  Theorem  3.12  is 

derived  from  condition  (3.24).  By 

* 

q.(')  :=  p!  <-)/2  +  c  /n ,  i  =  1 . n  , 

II 

"  ‘ 

r  ( • )  :  =  q  ( • )  -  \  q.  (•)  >  0  , 

....  .  *  .  .  i=l  •  • 

ql  :=  qi  *  r'  qi  :==  V  1  =  2 . n 

we  have  found  p  ,  qi ,  i  =  satisfying  (3.24).  a 

For  the  special  case  of  the  Poisson  equation 

I’  Pp")  =0,  i  =  q(-)  =  0,  that  is 

n 

Ay  +  f(-)  =  -  l  y  (■)  +  f(-)  in 

i=l  xixi 

y ( * )  -  q(-)  or.  3fi, 

the  stability  result  corresponding  to  Theorem  3.2  has  been  proved  in  Pereyra- 
Proskurowski-Widlund  (241.  We  can  not  include  this  special  case  in  our  general 
Theorem  3.8,  since  (3.27)  would  correspond  to  ct  =  0  in  (3.24).  The  proof  given 
above  breaks  down  for  ct  =  0. 


(3.29) 


F„y  :  = 
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4,  Stability  for  k  i  4  under  less  stringent  conditions. 

For  k  ^  4  a  much  easier  stability  proof  may  be  given.  If  P  is  a  i  '  .  permuta- 
T 

tion  matrix,  then  P  B  F  and  B  contain  the  same  diagonal  elements  only  in  a  changed 
order.  Further, the  elements  found  in  the  same  row  . olumn  with  will  be  found  in 

that  row  and  column  which  are  c".  i  *.*•:  ’  —ton  •  e.  '-r  ... 


Lemma  4.1:  Let  P  be  a  permutation  matrix.  Then  B  =  (b^  )  is  diagonal  dominant  if 

T  - 

and  only  if  P  B  P  =  (b^  )  is  diagonal  dominant  and 


min  { b .  .  -  Y  lb.  .  |  } 
1=1  j=l  J 


min  {b.  . 
i-i  11 


;  lb. .  | 

j-i  13 

j/i 


Further  we  need  the  elementary 


Lemma  4.2: 


Let  A  =  I 


3  =  1 


(a  ) .  let  each  nontrivial  row  in  B.  =  (b..  ,)  be  diaqonal 
i  £  -  1  It*  - - 


dominant  with  positive  diagonal  elements.  Further  let  for  each  row  of  A  exist  a  non¬ 
trivial  row  with  the  same  index  in  at  least  one  of  the  B ^  .  Then  A  is  diagonal  dominant 


and 


U  V  (  P  f  u 

min  {a.  .  -  T  la.  j  }  •  min  /  min  {  b.  . 

.  i=1.  ti  Jil  '.n  -  i=-v  •  31 

pyi  b. . . 5^0 

311 


y  ib,..!/ 


f=i 

t/i 


Now  if  we  want  to  prove  stability  it  is  enough  to  show  that  the  matrices  defined  in  (2.14) 
are  diagonal  dominant.  For  that  purpose  we  need  less  stringent  conditions  as  in  S3  (see 
(3.5))  . 

Lemma  4.3:  Let  in  (2.9)  a^,b^,c  satisfy 

a.  ,b. ,c  €  C  (fl) ,  0  <  a  <  a .  (• ) 

r  i  —  —  l 

(4.1)  )  o  <  Q*  <_  q(- )  :=  <=(•)/  \  a.(d  2Q 

i=l  1 

b. (0  „ 


a .  (•) 


<  P. 

—  l 


1 , . . .  ,n 


Then  for  h  <  min  {2/P^  >  and  k  =  1,2, 3, 4  the  matrix  B  =  )  in^(2.14)  is  diagonal 

dominant  with 


t  a  2 

(4.2)  min  {b.  >  |  b .  ,  I  }  >  c  h  ,c  >0. 

.  _  n  .  .  i]  1  —  +  + 

1=0  3=0 

j/i 
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Proof :  Under  the  assumption  (4.1)  and  for  h  min  'Z't 

n  i=l  ± 

c(-)/  I  a.  (•)  and  p  ( - )  :  =  b.  (*)/a.  (•)  (see  (2. ID) 

.  ,  l  li 

i=l 

2  +■  h2q  -  ( |-1  -  ~  p  /  +  ;  -1  +■  ~  p.  )  =  h“  q  ^  h“ 


Further,  the  — —  and  — 1  ,  v  =  have  alternating  sign  for 

*o  ?0 

—  >  0  .  So  Wv  :i:.i  ler  The  ars  ar.pt  1  -r.  i:.  :  :  :  < 


r:/  ■■  f  (1  +  j v  ‘  (|_1  + 1 Po +  ^  (1  +  r  V  +  •••  + 

k 

..  .  h  ..  h  .  ,  1  +  *2  “  ‘3  +  ‘4  '  lk 

-  +:-  j p0  +  (1  +  ipo'- - n - 


for  all  h  ••  0 


gk(s)  1  + 


3i  +  a2  ~  a3 


0  for  0  <  s  ■  1  . 


Now  one  verifies  either  by  straightforward  discussion  or  by  using  computers  that 


min 

0  <  s<  1 


gk(s)  >  0 
<  0 


for  k  =  1,2, 3, 4. 
for  k  =  5,6. 


So  (4.2)  holds  just  for  k  =  1,2, 3, 4  and  the  Lemma  is  proved. 

Combining  Lemmas  4.1  to  4.3  we  obtain  by  Gerschgorin '  s  Theorem  (21); 

n 

Theorem  4.4:  Let  (4.1)  be  satisfied  and  h  '  min  (2/pM  .  Then  the  matrix  A  is  roqula 

i=l  1 

and  there  are  positive  constants  Dk  such  that 
llA~1ll  <  D,  •  h"2  . 

In  contrast  to  Theorem  3.7  we  have  an  estimation  for  *'  -  - 

CO 

Since  in  the  proofs  of  our  Lemmas  4. 1-4. 3  the  coefficients  in  (2.1)  only  had  to  sati 

(4.1)  and  since  we  obtain  stability  for  the  nonlinear  problem  simply  by  proving  stability 
for  the  linear  problem  we  can  generalize  our  Theorem  ■• .  -  immediately  to  the  more  general 
case  given  in  (2.16),  (2.17): 
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Theorem  4. S ; 


(4.3) 


Let  2  be  the  unique  solution  of  (2.16),':. 


<  -■•  /  -  /  p* 


},  a:..:  le~ 


0  <  ?  1  ai<-,z(-)  ,zx  (•))  in  fl  ,  i  =  1,2,...  ,n  , 


o  <  Q*  1  {£<0'1'°---0)  (.,2(.),7z(-))  -  ‘i  a.t0'1'0) 


i=l 


(*  ,Z  (•  )  ,2  (•  )  )z  ( 

X  .  X  X 

1  11 


r  -1 

x(  l  a  ( •  ,z  ( • )  ,z  (•)))  <  Q*  in 

4  _  1  X  •  - 


i=l 


and  ,  with  P-j  in  ( 4 . 1 ) ,  let 


{  |{f(°i^i'°"-0,(-,z(.),Vz(.))  -  a(0'°'1)(.,z(-),z  (•)  )z  (•)■• 

'  i+2  1  x-  x.x. 

i  li 


X  (ai  (-,z  (•)  ,zx  ( * ) ) )  -1 1  <_  p?  in^  fl  , 
i  1 

with  continuous  functions  a.,  f  ‘  ‘  ‘ •  ,z  ( • )  ,  Vz  ( •) )  and 


.(o,...,qi,o,...,o) 


i+2 


(-,z(-),Vz(-)) 


Q  . 


Then  the  discretization  of  F  _in  (i.lG)  is  stable  in  a  nelsitfrourhood  o:  the  TOt 
solution  z. 
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5.  Convergence  and  asymptotic  expansion. 

With  the  stability,  proved  in  .3  and  4,  we  obtain  convergence  and  asymptotic  expan 
sions  by  studying  the  local  discretization  error.  Before  we  do  so  we  need  some  formal 

notations.  With  the  spaces  E  s.  C2(.:>  -  CL.)  and  E°  :=  CE)  *  C(  •.:)  and  the  grid 

>t.  •  '  •  •  .  •  •  I  .  »  s  ;  I  .  •  » ,  w  :ef 

■  pr:  ’ 


E. 


0 


h  :=  Eh  •  ' '  h  ,n 


.:)  (G, 


with  one  of  the  norms 

„0 


!ln  H 


Vh.- 


"Vh, 


X  £  D, 


h,n 


I  V  (x) 

h 


?  .)  •  »  I 


h,' 


(5.1) 


'V'h.2  :=  "  "h" h ,2 


[  h" 


I  :v*>^ 


1/2 


X  £  D, 


h,.: 


and  the  restriction  operators 


:=  f  '  '  ^  ,  »°  :  f  h 

l.  y  •'  y!n  11  5.  (u,v)  •• 


(u|„ 


h.» 


h, 


1  c.  n 

h,r. 


Now  we  have  (for  3  see  (2.1)) 

Lemma  5.1:  Let  in  (2.2)  the  solution  z  £  C‘ ),  — 


v.  v 

S:  1 . . ,3z  " 

n 


:.:f<0,C,V1 . Vn) 


‘1 


and 


v,  v 

3z.  X....3z  n 
1  n 


f(  •  ,  • ,  •  )€  C  (fi  x  B )  for  v.  +  . .  .+v.  =  v  =  1,.  and  let 


the  usual  formal  differential  operator  in  the  multivariate  Taylor  expans io:. 
n 

be 


given  as  (  \  4  ^~)v.  Then  the  local  discretization  error  <  f ) .  z ) 

i  =  l  “i 

is  given  in  regular  points  as 


^h‘(ip, 


n  q+1  2j  ...  j  ,  . 

r)(A  z)(x)  =  l  {-  a  (x)  [  —jr  - ) 

i  =  l  j  =  2  1 


v  ’5 


3x. 


(5.2)  { 


.u2  p  1  <  ?  ,  v  h2i  32j  +  1z(x).  3  1V  ,,  .  .  „  , 

+h  l  vi  \  /  1  r-rrT!  — — )  — 

V-1  v  1  i  =  l  j  =  l  -  ’  3xj  iJ 

+  0(h2(q+1)+a)  for  x  €  a  . 


qt  ? 

In  irregular  mesh  points  l  and  i  tsunmmg  the  derivatives  of  z,  have  to  be 

3=2  j=l 

changed  in  the  following  way: 

Let  z*  £  c2^q+*  ^+a((2* )  be  a  smooth  extension  of  z  ,  where  ft*  u  ft  is  large  enough 


to  contain  all  grid  neighbours  for  irregular  points.  Then  use 
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qvL  2h2  j2jz(x) 


(2j)!  2j  k+1 

jX 

i 


k+i  :.k+1z*(r.) 
h  1 


q  .2  j+2  2 j  + 1  s. 

Y  h  2  *(*>  i 

(2  j  +  l) !  ,x  2  j  +  1  +  2  (k+1) 


.  k+1  - 

s.  ,  _  -*  z  *  ( :, .  ) 

1  k+2  _ i_ 

k+1)  h  k+1 

‘x . 


with  £  , .  on  the  gridline  through  the  point  x 
in  e.  direction  and  i-th  co-ordinate  in 

—  i  — - - — -  — - - - - - 

(x  • (k-1) he . ,  x  +  he.) 
k.  11 

q+i  q 

for  the  corresponding  }'  resp.  V  in  (5.2). 

j=2  j=l  ~ 

Proof :  l-todiatviy  we  1  ..  )  wi  • 

u  (x+h)  -  u(x-h)  ...  :t1  h2-'  <2  j  +  1 )  ,  ,  2  .;-l+u, 

- - -  u  x  +  \  T2FTYT  u  (x)  t  0(h  ■  ) 

3  =  1 

and 

(5‘4)  \  t  ulx-t.)  $  2h2  U"X)  (21)  .  ,  _ 2 • : -2  +  . 

5 - — - -  =  u"  (x)  +  >  — — -  —  u  J  (x)  +  0  (h 

h2  j=2  (2])! 

for  any  u  ■  C2  ‘ (x-h ,x+h) . 

Now  z  e  C2(q+1,+a(S„  may  be  extended  to  a  z*  t  C2  *q+1'  +  a(2*) .  Such  extensions 
exist  ,  but  they  usually  do  not  satisfy  an  extended  differential  equation  (2.2). 

Further,  since  we  can  split  (5.2)  into  terms,  each  one  only  involving  derivatives  with 

respect  to  one  variable,  we  may  coniine  the  iiscusr  Ion  now  to  -.he  variable- 
Since  we  need  trie  iocal  !  Inc  re- in  +  1  1  m  error,  wo  have  :■  r,;  lace  y.  in  by  =  .  , 

ot  tained  as  ?.  j=P^(x+he.  ),  where  i_  is  define  i  P,  (  >:-(  v-.  Ib.e.  )  =  s,  ,  v*i,...,h, 

Pj,  (x* /  -  g(n.  )  .  .0  we  nave,  :  •!.<•  well  know:.  error  formulas  ‘or  polynomial 

interpolation, 


(x+he .  -x*)  (x+he  .  -x) .  . .  (x+he  .  -x+  (k-1)  he  . )  „  , , 

z  (x+he  . )  -  *  =  - - i - i - i-  •  z  (k+11(U 

(k+1)  ! 

,k+l 

Si  .k+1  3 

k+1  3x . k+1  ' 


By  inserting  (5.5)  into  (5.4)  we  evidently  obtain  (5.2),  (5.3). 


-  ( p4  - 1 ) 

Combining  formulas  (5.2)  and  (5.3)  and  observing  that  there  are  only  0(h  ) 

irregular  points  one  obtains 

Lemma  5  ■  2:  The  local  discretization  error  (^>^r)  ( A^z )  admits  an  asymptotic  expansi : 


the  form 


lx 


for  2q  <  k-1  +  1/a,  and  0=2,  k  <  sr  c  =  ",  -  :  •  , 
whe re  the  g^(  are  defined  in  (5.2),  (5.3). and  q  <  q  as  in  ^enrria  5  ■  :  ■ 

Again,  since  (5.1)  can  be  split  into  the  terms  mentioned  in  (2.11),  we  may  use  t!.< 
arguments  in  [7]  to  prove  Theorem  5.3.  In  this  Theorem  we  only  give  conditions  for  a  ,  f 
and  g,30  to  ensure  z  e  C2q+2+a(li).  The  more  or  less  "geometric"  type  conditions,  given 
e.g. ,  in  Agmon-Douglis-Nirenberg  [1],  §7  are  only  referred  to.  If  they  are  valid,  the 
smoothness  results  for  the  e^ ;  are  ensured. 

Theorem  5 .3:  Let ,  for  k  <4  and  k  <_  6,  the  II-  11^ -norm,  a  =  °°  and_  0=2  »  the  condition 

(4.3)  and  (3.26  )-(  3.27)  or  (3.58)  be  valid  respectively.  Let  further  ;:£  C 1  ( H ) ,  a  -  £  1  'V.i, 

f  (-,-,-)  c  C2q+a(!)xB)  ,  30  €  C2q+2  +  a,  g  t  C2q+2+a(30)  for  0  <  a  <  1  and 
k  -  1  +  l/o  2q+a,  and  "geometric"  conditions  (see  tl]  ,  §7)  be  satisfied. 

Then  z  €  C2q+2+a(fl)  and  we  have  for  the  approximate  solution  5 


k-l+l/a 


5h  -  V  -  fih  J  h  e 21  +  0(h 


for  q  =  max{m  t  3N|2m  <  k-l+l/o,  m  <_  q)  . 

here  0(hk  ^+^°)  refers  to  the  difference  with  respect  to  II  -II  t0=2>°°-rurt:-.<n  we  nav 
e2  t  ( *)  e  C2(q  ')+2+a  (0)  ,  e2  ^  (• )  =  0  on_  3fi  ,  1  =  1,2  ....  ,q  . 


(5.7)  corresponds  for  (3.29)  to  the  result  of  Pereyra-Proskurowski-Widlund  124]. 

A  result, given  in  Agmon-Douglis-Nirenberg  [1],  §7,  may  sometimes  help  to  avoid 
difficulties  arising  from  the  violation  of  30  e  Q2<3+2+a  ^n  isolated  points.  Since  in 
all  our  proofs  we  only  need  differentiability  properties  along  grid  lines,  one  only  has 
to  choose  the  grid  such  that  its  lines  avoid  non-smooth  boundary  points.  Then  the 


results  in  [1],  §7  may  be  used  to  obtain  results  for  the  case  hQ  >_  h  >_  h^>  0  ,  which 

might  sometimes  be  useful. 


I 


6.  Two  special  cases. 

The  breakdown  of  the  asymptotic  expansion  in  Theorem  5.4  is,  for  a  sufficiently  snoot! 
situation,  caused  by  the  interpolation  error,  given  in  (5.5).  This  error  will  be  less 

disturbing  in  one  of  the  two  following  cases: 

(6.1)  Let  s.  =  0,  i  =  p  (see  (2.6)),  in  (5.5)  for  all  x  c  si  .  . 

l  1  m  h,i 

(6.2)  We  admit  only  h,  such  that 

s.  =  const,  (independent  of  h  ),  i  =  u . ,  p  ,  in  (5.5)  for  all  x 

i  1  m  h,i 

The  second  case  is,  admittedly,  somewhat  artificial,  but  if  it  is  satisfied  we  have 

instead  of  (5.7)  the  relation 

C  -  A  z  =  A  J  h2le  +  0(hk  +  1/°) 
h  h  h  ' .  2 1 

1=1 

(6.3)  /  ^  (6.2)  is  satisfied,  with  k+l/o  <  2q+a  and 

q+  =  max  {me  jn  |2m  <  k+l/a,  m  <_  q  ). 

In  some  cases,  when  (6.2)  is  satisfied  we  may,  by  a  proper  choice  of  (different)  stepsizes 
in  the  directions  of  different  co-ordinate  axises.even  obtain  si  =  0,  that  is  (6.1).  This 
case  has  been  treated  in  Pereyra  [23]  for  uniformly  elliptic  operators  of  the  following 
(casi linear)  type 


(6.4)  F  y: 
c 


l  *  (#)  +  f(*.y(*),7y(*))  in  15 

i,j=l  3  *i  j 

y ( * >  -  g(-)  on  3(5 


r  2 

with  l  a.  .  (•)£.£.  >_  a  II  £ll  in  (5  with  a  >  0  and  for  any  real  vector 
i ,  j=l  13  1  3 

T 

£  =  (£  ,...,£  ).  The  discretization  of  (6.4)  under  the  condition  (6.1)  is  obtained  by 

in  J 

using  (2.4)  in  regular  and  irregular  mesh  points  and  by  replacing  in  irregular  mesh  points 
the  n^Cxihe^)  with  x±h  ei  e  3(5  by  the  boundary  values  g(x±he^).  Then  Pereyra 
1231  has  proved  (the  smoothness  properties  of  the  ^  in  (6.5)  may  be  checked  by  some 
additional  computations) . 

Theorem  6.1:  Let  (6.1)  be  satisfied  and  the  operator  F^  in_  (6.4)  be  discretized  as 
described  above.  Further  let  a_  e  C2q+“(r<),  f  (•,•,•)  e  C2q+a  (j)xB)  and  3(5  e  c2q+2+a 
g  e  C2q+2+a(3(5)>  z£C'*'(S5)  and  "geometric"  conditions  for  3(5  be  satisfied  .Then  z£  C2q+^a(fi)ar 
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7. _ D iscre t e  Newton 


and  deferred  correct  ion:; 


As  an  immediate  consequence  of  the  results  in  w»*  m.n  it. •»*  .u.i  on  •  si  «  i 

polation  techniques  to  get  higher  order  results.  The  general  .u-ptwi.l;  i  .iiv«*u  m 
Bulirsch-Stoer  1131  and  applied  to  the  Tiriehlet  problem  for  t:v  roi;;::on  equation  m 
Pereyra-Proskurowski-Widlund  (241.  So  we  only  give  the  result  here. 

Theorem  7.1:  By  applying  Richardson  extrapolation  once  resp.  twice  to  the  methods  given 
in  Theorem  5.3  for  k  *  4  resp .  k  =  6  we  get  approximations  of  order  3.5  resp.  5.5  in 
S  -norm*  If  (6.2)  is  valid  orders  4.5  resp .  6.5  are  obtainable  by  •  :•  „  1/ 

three  extrapolation  steps.  Even  higher  orders  (up  to  2q+l,  see  Theorem  6.1)  are  obtainable , 
if  (6.1)  is  satisfied. 

Using  Richardson  extrapolation  requires  the  solution  of  the  nonlinear  system 
(*hF);,  -  0  for  decreasing  h  .  That  implies,  for  higher  dimensions  n  %s  *  >.  v :  v .  .  ‘  •  I 

sizes  c_  the  sys terns  or  nor.l  Ir.esi*  er  iar i  'r.r  tr  oe  snivel.  aie  able  to  avtl  *  .  ’ 

disadvantage  by  either  applying  Pereyra’s  (22]  iterated  deferred  corrections  or  iterated 
defect  corrections  or  descrete  Newton  methods  (6,8)  . 


The  first  two  methods  require  for  every  correction  step  the  solution  of  a  nonlinear 
equation  (of  the  same  size)  where  the  original  right  hand  side  D  has  to  be  changed.  Ir. 
discrete  Newton  methods  only  the  original  nonlinear  equation  has  to  be  solved. 

The  corrections  are  obtained  by  solving  linear  systems  (of  the  same  size)with  a 

fixed  matrix  and  different  right  hand  sides  and  therefore  need  much  less  computational 
time. 


Most  of  the  conditions  for  the  applicability  of  iterated  defect  corrections 
are  easily  verified,  whereas  for  discrete  Newton  methods  we  need  a  lot  of  technical 
formalism.  So  we  prove  the  results  for  iterated  defect  corrections  and  present 
discrete  Newton  methods  and  sketch  deferred  corrections  without  proof.  For  deferred 
corrections  we  have  to  approximate  (only  k=6  is  discussed  here)  the  local  discreti¬ 
zation  error  given  in  (5.2)  as  (gyXA^HxXh^tx)  +  h4b.,(x)  +  ... 

Now  we  have  from  (5.2),  with  — — — - .  f:=f^  ’  ,vl*”’vn^ 


3z^'ll. 


.  3z  n 
n 


••+vv. 


0  v  »* 

h  bl(x)  ={-£{■ 
i  =  l 


ai(x)  34z(x)  33z(x)  3  r 

°  +  - ^92“  ^x»2* x),7z(x))} 

i 


3x.‘ 

1 


3x. v 
1 


and 
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h4b2(x) 


v,1*  |  n  a. 

=  120  1  Jx  {‘ 


a. (x)  .6  .  ,  5 

3  2 ( x )  3  z(  x )  3 

6  ~ — 5 —  f(x,z(x),V  z(x) )} 

9x.  3x.  i 

i  i 


+  3  (  .1  ~  Z3X)  317 }  f(x,z(x),V  z(x))l 

1“1  ox.  l  J 


2  2  4 

and  so  we  need  approximations  for  h  (x)  an;  h  (x)  +  h  t>2  (x)  of  the  order  4  a 

6  respectively.  Now  the  same  remark  ist  appropriate  for  irregular  points  with 
respect  to  s  4  1  for  the  behaviour  of  the  ct^/a^  as  given  in  Fereyra-Proskurowsk i- 
Widlund  [24],  p.  9. 

In  iterated  defect  corrections  and  in  discrete  Newton  methods  we  have  to 
discreti2e  equations  of  the  form 

(7.1)  F^y:  =F  y+d,  with  d:  =(%in  . 

d  oJ  \e*on  dUJ 

For  defect  corrections  F  is  given  as  F  in  (2.1),  whereas  F  £  l(C2(ft)  n  C(£2 ). 

o  o 

C(ft)xCOS)))  for  discrete  Newton  methods,  with  the  discretization,  defined  in  2.5 
we  obtain 

'  (W  nh  a(Vo)nh  +  V’ 

(7.2)  ■  where  ip,F  is  described  in  §2  and 

h  o 

V  L(C(Q)x  C(3fi),(rh  nnSJ)x(Gh  J13K)). 

For  the  considerations  below  we  have 

(7l3)  d  Z  (o  on  3fl  )  and  therefore  V  =  (^b.n™)’ 
whereas  for  e**0  in  (7.1)  the  component  of  f^d  in  ^0(2  is  more  compli¬ 
cated  since  the  boundary  values  enter  the  difference  equations  for 
irregular  points  in  (2.6).  However,  in  either  case  S^d  is  independend  of  t^. 
Discretizations,  for  which  (7.1)  implies  (7.2),  are  called  strongly  linear,  see 
Bohmer  [8].  In  [8]  we  have  proved  that  ,  if  Fy  t  d  =  0  is  uniquely  solvable  for 
small  enough  ]|  d  ||  and  if  the  strongly  additive  discretization  method  generates 
(l^d  independend  of  ri^,  a  stable  <p^F  and  an  error  asymptotic,  where  the  coefficients 
e2i  dn  are  elements  of  linear  spaces  with  certain  properties,  see 
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Method  1:  Continuous  piecewise  polynomials: 

With  a  *  >  k,  k  as  in  (2.5),  define  7  as  (TT^  are  the  polyr.o 
1  and  Is]  -  entire  part  of  s) 


(  P)  *  “  iTV  I  /  y 

x  (xuX,x(u+i)>. 


(7.5) 


e  nx,u=o,i,...,y*:  =  [  y ]-  -, 


P  (>:.)  =  n.,  i  =  yX, . . . ,(y+i)X, 

y  i  i 


P  *  Try  /  \ 

y  +i  h  x(p%i  )X,xrtl 


P  .  (x.)  =  n.,  i  -  r+l-X, . . . ,r+l.  a 

Li  +1  1  1 


Method  2:  Spline-functions: 

Define  7  as  a  spline-function  of  degree  X  >  k  as  discussed  ir.  [4,6,8] 


Method  3:  Overlapping  polynomials  (discrete  splines): 
Let  X=X*+2  K  >k,  <  =  1,2, 3, 4  and  define  T  as 


e  il 


(  P_:  '  Trlv,  I  *-  “X 

°  h|(x  x  )  A> 

xo,XX*+<  ’ 


(7.6) 


P0(^v)  =nv>V  :  0,...,x, 


P  :  =  Try  - 

V1  hI(x  *  x  *  )  * 

^XyA  +K,X( y+1 )A  +k' 


e  n.,y  r  i . ££]-  i. 


P,  (x  )  =  n„  ,  V  =  yX*, . .  .  ,(y+l)X*+2ic,  y  =  -1. 


y  v 


\  The  last  polynomial  is  defined  analogously  to  (7.5)  □ 
Method  4:  Symmetric  formulas  and  polynomial  extrapolation  [  8,24]  : 


Interpolate  the  points  v  =  0,l,...,k,v  =  r-ktl , . . . ,r+l , 


v  =  1,...,X  t  1,  and  v-  r-X,...,r  by  polynomials  P,  .  ,  P,  ,  P,  .  and  P, 

Kji  K*e  A 9 1  a  ,e 

respectively,  and  compute  extrapolated  provisional  values 


(7,7)  %  =  Pk,i(xv),V  =  and  %  z  Pkje(xv)’  V  =  r+2’r+3 . 


and 


(7.8) 


nv  =  PA,i(V’v  =  O’-1*-2*”-  and  V  PA,e(xv)>v  =  r+1’r+2 . 
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tc  u  ) 


Ter  irregular  points  x  6  .  the  changes  corresponding 

*•  »  - 

mate .  in  the  z  (x)  s  r..  ( x )  are  known,  the  V.t  (x>  may  he  comcute;  hy 

o  r. 

or.e  of  the  methods  iefir.ed  above.  In  [i , r]  we  nave  sr.owr,  tr.at,  unoer  c-or.  uti 

satisfied  here,  we  do  not  need  ir.  f*(z  )  exactlv,  but  that  we  rr.av  nn: 

over  to  an  f*(  )  which  is  obtained  fror.  F'(z  )  by  replacing  Vz  ( x  >  by 

z  (x+ne.)  -  z  (x-he.)  z  (x-t-he  )  -  z  (x-he  ) 

(7.:.?)  i— - - - r-2 - —  —  - —  )»Vzo(x). 


how  we  are  reahy  tc  formulate  our  discrete  howto:,  tat  he  :  for  the 

probier,  I,..  ).  Let  iP.,F,  at  refine.;  in  (.'.*•  .  F  .  *(  z  > : -T'  ( z^  )  a:.; 

Fj*(z  ):=  r*(z  ),  by  replacing  ir,  F'(;a)  V  -  <  >  .  :  :  r 

(7.13),  (ft  ff  ( :•  )  as  in  (7.11)  and  use  one  of  the  operator.'  1  in  "etho:  i — . 
n  v  o 

/  (i)  giver:  the  starting  value  r,}j  q  such  that 


(7.14) 


((p  r )C.  =  i, 

n  n,o 


t  c  ,  :  =  Tr  ana 

t-l  h,».-l 


compute  ;h^,e  =  1,2, . .  :  - 


n 


(ii)  (cnr»(z^  ))(r  -r  „  )  =  -  (  -  \ 

h  °  h’*  ■1’^1  10  on  G,  n  80 ' 

h ,  n 


This  method  is  called  discrete  Newton  method.  Ac  indicated  already  in 
Pereyra-ProskurowsFi-Widlund  [24]  the  order  results  for  deferred  corrections, 
als  well  as  for  discrete  Newton  methods  and  iterated  defect  corrections,  are 

K  “  1.  ^  * 

rather  poor  due  to  the  fact  that  in  (5.7)  we  have  only  0(h  r-  ), 
k+1 

instead  of  an  0(h  )  which  would  correspond  to  the  well  known  order  results 

k  =  0  and  k  =  1 . 


Theorem  7.2:  Let  k  =  6,  >■  >  G  in  Method  1-4,  and  let  the  conditions  in 

^  3+cx  - 

Theorem  5.3  be  satisfied  with  q  -  3  and  a  >  0.5,  so^  z  €  C  (H).  Further, 
let  be  defined  either  as  the  first  deferred  correction  or  by  (7.10)  or 

(7.14) ;  Then 

(7.15)  C  =  A  z  +  0(h3,5). 

h,l  h 

Further,  if  (6.1)  is  satisfied  we  have  instead  of  (7.15) 

(7.16)  C  ,=  A,  {z+hV  .}  +  0(h4,5) 
h, 1  n 
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?or 


If  we  compare  the  different  types  of  defect  ar.d  feferro.:  correct  I 
with  the  Richardson  extrapolation,  the  main  advantage  cf  the  extract la'i 
to  work  up  to  order  5.5  compared  to  order  3.5  for  the  orther  methods  an:  c 
provide  us  with  asymptotic  upper  and  lower  bounds  for  the  exact  solution 
(see  Bulirsch-Stoer  [13]  ).  However,  one  has  to  pay  for  these  advantages 

by  much  higher  computational  time  than  for  discrete  Newton  methods  or  one 
of  its  equivalents.  Expecially> for  nonlinear  problems  (2.1)  the  discrete 
Newton  methods  yield  the  corrections  in  a  relatively  small  amount  cf 
additional  computations,  by  solving  linear  equations  with  a  known  matrix  t 
(see  8.).  For  those  problems  the  discrete  Newtor.  methods  are  ever,  superior 
iterated  defect  and  deferred  corrections  which  furnish  the  corrections  : v 
solving  nonlinear  systems  of  the  initial  size  with  very  good  known 
starting  values,  whereas  Richardson  extrapolation  needs  the  solution  of 
nonlinear  systems,  with  numbers  of  equations  proportional  to  1/h  '. 


8.  Computational  remarks: 

In  equation  (2.6)  the  fraction  j/a0  ^  will  become  verv  large  if  s , 
is  close  to  1,  that  is,  if  the  irregular  point  x  €  .  is  close  to  the 

boundary.  Therefore  we  will  scale  the  equations  in  (2.6)  always  so,  that 
the  coefficient  in  the  linear  part  of  n^(x)  becomes  2na. 

If  (2.2)  is  linear  we  may  use  any  good  method  to  solve  the  linear 
system  (2.4)-(2.6).  The  choice  of  the  method  very  much  depends  on  the 
properties  of  f,  see  e.g.  Rosser  [25]  .  In  many  cases  it  will  be  possible 
to  get  the  corrections  in  the  discrete  Newton  method  (we  choose  v=  2  in 
(7.10)  for  such  a  case)  by  considerably  less  computational  effort.  That  is 
for  example,  true  if  the  differential  equation  in  (2.2)  is  the  Poisson 
equation  (treated  in  Pereyra-Proskurowski-Widlung  [24]  or  the  Helmholtz 
equation  Az  +  cz  =  0,  c  =  constant) . For  both  equations  one  can  very  well 
use  the  fast  Laplace  solvers  as  described  in  [24]  and  then  gets  the 
corrections  pretty  cheaply. 
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If  (2.2)  is  nonlinear,  the  system  (2.4)-(2.S)  will  be  r.cr.lir.ear  too 
most  cases  one  will  then  use  the  Newton  method 

(8.1)  V,,lVi,(VV1)*-V,V1'  *=1,2,... 

to  obtain  the  solution  £  *=  lim  n  of  (2.4)-(2.8).  If  an  appropriate  stoping  criterion 

"  £■+«>  £ 

for  the  iteration  (8.1)  has  been  chosen,  so  C.  :■  r\n  ,  one  may  often  use  the  alreadv 

'  h  Jc 

known  matrix  (tf^F ) 1  (n^  ^)  for  (F*(z  ))  ^  (7.10).  So, again, the  corrections  by  the 

discrete  Newton  method  are  obtained  relatively  cheap. 

Finally,  we  have  pointed  out  already  in  7.  that  for  a  proper  choice  of  h  the 

discrete  Newton  method  may  be  executed  on  a  fixed  grid  r.  and  (6.2)  will  usually  be 

h  ,n 

satisfied.  So  the  chances  are  pretty  good  that  a  r -rrec'lor.  improve.  • 

mation.  In  that  case  it  is  worth-while  to  use  A  =  k+1  in  Methods  1-4. 

If  (6.1)  is  satisfied  ,and  the  situation  is  smooth  enough  ^it  is  advisable  to  take 
pretty  big  (X  >_  10),  Because  of  the  loss  of  asymptotic  terms  in  the  methods  discussed 

in  Theorem  7.2  Richardson  extrapolation  might  still  he  ir.te:  i:.g  for  lest 

smooth  situations. 
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j.  Numerical  examples 

As  demonstrating  examples  we  use  boundary  value  problems,  all  the  ?a-e 
domain  fl,  defined  by  the  differential  equations  F.z  =  0,  i  =  and  bv 

the  boundary  values,  which  are  the  restrictions  of  the  exact  solutions 
z.,  i  =  1,...,4,  to  the  boundary  Si;.  So  we  have 

f  R:  =  { ( x ,y  ) !  ( x-0. 5  )'  +  ly-1)2  <(0.36)*}  ; 

r.z:  -  z  +z  +  Zz  -  1,  z(x,y)  =  sin(xy); 
u  xx  yy 

,  r,  ,  ,  F.z  =  z  +z  +  13z  =  C,  z(x,y)-  sindx+tv); 

( 9 . 1 ) X  z  xx  yy 

F,z  =  z  +z  t  ie  zZ=0,  z(x,  •)*  lr.(x+v); 

3  xx  yy 

'l  F,  z  —  z  tz  -lix‘=0  ,  z<x,y)  -  (x+y)  *. 

\  4  xx  yy  ^ 

With  the  discretization  in  §2  we  obtain  b  F.  and  b,  F,  as  systems  of  linear, 

n  1  n  c 

bhF3  and  b^F^  as  systems  of  nonlinear  equations  respectively.  The  first  step  it. 
(7.10)  and  (7.14)  is  to  solve  (b^F. )C^  =  0.  The  nonlinear  systems,  for  i  -  3,-, 
are  solved  via  the  usual  Newton  method, see  it, (which  is  not  to  be  mixed  up  with  the 
discrete  Newton  method  below).  The  matrices  of  all  the  problems  above  are 
positive  definite  and,  for  the  special  case  of  Ct  in  (9.1),  symmetric.  So  we 
solve  the  linear  equations  by  using  SOR  methods  combined  with  techniques  to 
estimate  the  optimal  relaxation  parameters.  For  the  linear  problems  F^  and  F, 
these  relaxation  parameters  have  to  be  estimated  only  once,  since  the 
discretization  b^F^,  i=l,2,  and  the  discrete  Newton  method  are  based  on  the 
same  matrix.  For  the  nonlinear  problems  the  matrix  in  the  discrete  Newton 
method  is  often  replaced  by  the  last  matrix  in  Newton' s  method  to  solve  the 
nonlinear  system  (bj^FJ^O,  i=3,4.  So  the  estimation  for  the  optimal  relaxation 
parameter  in  the  last  Newton  step  again  may  be  used  in  the  discrete  Newton  method. 
Therefore  the  iterations  in  the  discrete  Newton  method  are  rather  cheeply 
available . 

In  the  following  Table  we  give  the  errors  ]l  Ch  ^-^z  ||^for  v  =  0,1,2 
and  for  k  -  3,. ..,6,  where  k  is  given  in  (2.5  ).  The  defects  in  the  discrete 


■BPMPLW 
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Newton  methods  are  compute:;  with  methods  3  and  see 
use  overlapping  polynomials  with  <  -  2  and  ir.  method 
results,  where  we  have  usee  polynomials  of  degree  3  i 
provisional  values  for  the  outer  points  which  we  need 
divide;  differences.  Because  or  ('.lb)  we  have 


Ch,v-1  -  “r 


h,v-l  ’h,v 


Therefore  the  correction  v  .-5v  v  tiay  be  used  as  an  estimation  for  the 
error  Ch  especially  Ch  C-n  ,  has  been  user  to  estimate 

The  uunihers  q  are  defined  as  quotients  between  these  estimated  errors  and 
real  errors. 

The  following  numbers  are  computed  on  a  UN'IVAC  1108  with  double  prec 
(about  18  digits).  I  want. to  thank  cand.  math.  H.  Offerman.n,  who  did  the 
computations . 
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Example 

1 

Exampie 

1 

method  3 

method 

4 

method 

method 

4 

k 

V 

error 

q 

error 

q 

error 

q 

error 

q 

3 

0 

1.3, -06 

1.0 

1.3, -06 

1.0 

3. 2, -05 

1 .0 

8. 2, -06 

1.0 

1 

1.6, -09 

1.1 

2. 3, -08 

1.0 

3. 0,-07 

1.0 

4 .3,-07 

1 .0 

2 

1.7, -10 

0.6 

7. 6, -10 

0.8 

1.3, -08 

o.e 

4. 2, -08 

0.9 

4 

1 

4. 2, -10 

1.0 

2. 0,-08 

1.0 

1.4, -07 

1.0 

4. 3,-0? 

1 .0 

2 

8. 6, -12 

0.9 

2. 4, -10 

1.0 

2. 3, -09 

0.8 

3. 4, -09 

1 .3 

5 

1 

4. 3, -10 

1.0 

2. 0,-08 

1.0 

1.4, -07 

1.0 

5,3, -08 

1.0 

2 

1.2, -12 

1.0 

2. 4, -10 

1.0 

4. 6, -10 

1.0 

6  ,4  ,  -10 

2.0 

6 

1 

4. 2, -10 

1.0 

2. 0,-08 

1.0 

1.5, -07 

l.C 

6 .5,-08 

l.C 

2 

9. 6, -13 

1.0 

2. 4, -10 

1.0 

2. 8, -10 

1.0 

3. 8, -1C 

2.4 

Example 

3 

Example 

4 

3 

0 

2.5, -06 

1.0 

2. 5, -06 

1.0 

2. 4, -05 

l.C 

2. 4, -05 

l.C 

1 

6. 6, -09 

1.1 

2. 3, -08 

0.9 

7. 7, -08 

. 

3.2, -07 

1.0 

2 

1.1, -09 

0.6 

2. 5, -09 

0.8 

1.6, -08 

0.6 

3. 5, -08 

0.7 

4 

1 

2. 8, -09 

1.0 

1.8, -09 

1.0 

6. 4, -08 

1.0 

3.6, -08 

1.1 

2 

1.2, -10 

0.8 

1.6,- 10 

0.8 

2. 4, -09 

0.7 

3. 0,-09 

0.6 

5 

1 

3. 2, -09 

1.0 

3. 5, -10 

1.0 

7. 6, -08 

1.0 

4. 2, -09 

0.9 

2 

1.3, -11 

0.8 

2. 0,-11 

0.4 

2. 6, -10 

0.7 

1.1, -09 

0.6 

6 

1 

3. 3, -09 

1.0 

1.9, -10 

0.9 

7. 8, -08 

1.0 

3. 4, -09 

1.1 

2 

5. 2, -12 

1.0 

3. 4, -11 

0.7 

3. 2, -10 

1.1 

1.5, -09 

1.1 

Table 


From  this  table  we  see  that,  even  for  k  =  3,  the  first  correction  is  worthwhile 
and  for  larger  values  of  k,  especially  k  -  6,  two  corrections  considerably 
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